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^\ , The degree of ionization of a nondegenerate two-dimensional electron-hole plasma is calculated 

0^ ' using the modified law of mass action, which takes into account all bound and unbound states in a 

0^ I screened Coulomb potential. Application of the variable phase method to this potential allows us 

to treat scattering and bound states on the same footing. Inclusion of the scattering states leads 
f~| ' to a strong deviation from the standard law of mass action. A qualitative difference between mid- 

^ , and wide-gap semiconductors is demonstrated. For wide-gap semiconductors at room temperature, 

when the bare exciton binding energy is of the order of ksT, the equilibrium consists of an almost 
r~^ ' equal mixture of correlated electron-hole pairs and uncorrelated free carriers. 
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I. INTRODUCTION 



The drive for ever higher stajjage capacity has led to the development of semiconductor lasers operating in the blue 

. . spectral region, based on ZnSdd and GaNE Along with the large energy gap of these materials comes a lame exciton 

binding energy, of the same order as fcsT at room temperature. As has been well known since the seventies,a excitonic 

gain processes are important in wide-gap semiconductors, and their importance is further enhanced in quantum well 

structures where the binding energy may be considerably larger than ksT (e.g., ~35 meV in ZnCdSe/ZnSe quantum 

'^ ] wells). 

^ ' Theoretical treatments of GaAs and InP based lasers are well established using a microscopic many-body approach 

O I based on linear response theory.cl Screening and band-gap renormalization effects are included assuming the injected 

carriers form a completely ionized electron-hole plasma. Such treatments have been successful in explaining many 

of the observed features of mid-infrared laser diodes. Complex valence band effects and strain effects as well as 

—< , carrier thermalization effects have all been included at various levels of complexity. In this way, a relatively complete 

^ ' understanding exists for the basic operation of these lasers. 

In wide-gap semiconductors however, the strong Coulomb interaction leads to the existence of bound exciton states, 
which persist even at elevated densities and temperatures. As such, the conventional assumption that the inversion 
vQ I is in the form of an electron-hole plasma with no excitons present deserves closer examination. A self-consistent 
^^ I description where both bound and unbound states are treated on an equal footing is required. Unfortunately, as far 
^\ ' as we are aware no comprehensive theoretical treatment of this problem exists. Treatments based on bosonic exciton 
On 1 operators have been proposedd but this approach breaks down at high injection when the screening of the Coulomb 
potential weakens the binding and produces a population of unbound scattering states, which clearly do not exhibit 
bosonic character. On the other hand, a treatment based around fermionic electron and hole operators is complex 
when higher order excitonic correlations are important. Q'Ll A natural concept in considering this issue is the degree 
of ionization in the interacting electron-hole plasma, and in this paper we calculate this for a two-dimensional (2D) 
plasma. We focus on 2D for two reasons. Firstly, most modern semiconductor lasers are fabricated in quantum 
Q ' well heterostructures. Secondly, the presence of at least one bound state in the attractive 2D potential requires a 
^ nonperturbative treatment of the screened Coulomb interaction. 

sj ' We will be mostly interested in the plasma properties induced by the pair Coulomb interaction between charged 

• i-H . particles, neglecting band gap renormalizaUop and phase-space filling effects, which have been extensively studied in 
^^ ' both three-dimensional (3D) and 2D cases.ErQ These effects can be neglected only in the low-density (nondegenerate) 
JH ' limit, which is defined in 2D by the inequality 

nXlj/g < 1 , (1) 

where n is the 2D carrier density, g is the spin degeneracy factor of 2D particles, and Am = (27r?i /MksTY^'^ is the 
thermal wavelength. For the two-component plasma the lighter carrier (usually electron) effective mass must be used to 
evaluate the thermal wavelength to ensure that condition (Q) is valid for both types of carriers. Inequality (H) provides 
that the motion of excitons can also be considered as classical. For GaAs at room temperature Xm^ ~ 1.66x10"^ cm, 
electron spin degeneracy ge — 2 and condition (n^) is satisfied for n < 7.2 x 10^^ cm~^. The electron effective mass 
in wide-gap semiconductors is usually at least two times larger than in mid-gap semiconductors, therefore condition 



(||) is valid over a wider range of carrier densities (e.g. for ZnSe at room temperature inequality (||) is satisfied for 
n ^ 1.7 X 10^^ cm~^). Thus, the nondegenerate (Boltzmann) limit is not only a convenient approximation, in which 
the Coulomb interaction is not hidden by the band-filling effects, but it also gives a realistic picture of the electron-hole 
plasma in wide-gap semiconductors at room temperature and moderate carrier densitijes. Lasing at anomalously low 
densities (below the Mott density) has been reported in ZnCdSe/ZnSe quantum wells.ll5l 

Following anjapproach applied in 3D to nuclear matter Jiil an ionic plasmajlj and the electron-hole system in excited 
semiconductorsEJ we divide the total electron (hole) density between two terms: 

ua = n° + nr"^ . (2) 

The first term n1 is the density of uncorrelated quasiparticles with renormalized energies. This term is that part 
of the total density which is independent of the inter-particle interaction (see Appendix g). All correlation effects 
both in the bound and continuum states are incorporated into the second term n™'"' , which is called the correlated 
density. The lower index in Eq. (0) is a species index, a = e for electrons and a = h for holes. It is also useful to 
introduce the degree of ionization of the electron-hole plasma 



n? 



(3) 

which characterizes the deviation of the thermodynamic properties of the electron-hole plasma from those of the ideal 
gas (a = 1). The knowledge of the degree of ionization is essential in determining the dominant lasing mechanism. 
When a is close to unity the main lasing mechanism is stimulated-emission from the free-carrier plasma, for lower 
values of a several excitonic gain processes have to be considered.til 

In the nondegenerate limit there is no need to go beyond two-particle correlations. This allows us to separate clearly 
the role of the inter-particle Coulomb interaction from the phase-space filling effects. In this limit, the correlated and 
uncorrelated densities are related by 



E"^° 



—^— Zab , (4) 



where (3 = l/(fc_Br), /Xafc = MaMb/{Ma + Mb) is the reduced effective mass, and Zab is the two-body interaction 
part of the partition function. This relationship is derived in Appendix g. Note that due to charge-neutrality the 
total electron-hole density Ue = rih = n is independent of species, whereas n^ ^ ri° and n™'"'' ^ n™'"'' if electrons 
and holes have different effective masses. 

The electron-hole part of the partition function which exhibits bound states (excitons) is given by 
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where mh is the projection of the angular momentum onto the axis normal to the plane of 2D motion (m = 
0, ±1, ±2, . . .), h k^ /{2iJ,eh) is the energy of the relative motion of the unbound (scattered) elecimn and hole, hk 
is the magnitude of the relative motion momentum, (5m (fc) are the 2D scattering phase shifts JlJH Em.u are the 
bound-state energies (index v enumerates bound states with given m), and the double Siifn in the first term ranges 
only over bound states. Equation (^ is the 2D analogue of the Beth-Uhlenbeck formulaO and it can be derived in 
the same fashion as in the 3D case, as shown in Appendix H. 

The scattering (integral) term in the right-hand side of Eq. (^ gives the contribution to Zeh of the continuum part 
of the eaetra^ spectrum. This term is usually neglected in calculations of the ionization degree of the electron- hole 
plasmaH'tirta In what follows we will show that at high enough temperature the scattering term is comparable to 
the bound-state sum and indeed this term has to be taken into account to ensure continuity of the partition function 
whenever bound states disappear with increasing screening.l23 

The electron-electron and hole-hole parts of the partition function Z^e and Z^h contain the scattering term only. 
To calculate Zaa one must take into account the Pauli exclusion principle for identical particles, which modifies the 
sum over m. The electron-electron (hole-hole) part of the partition function is given by (see Appendix |B|) 



Zaa = I yi2- (-1)-} r'-^e.^ (-P^) dk . (6) 
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Here we assume that both electron and hole states in quantum wells are two-fold degenerate. The only difference 
between Z^h and Zee arises from the difference between electron and hole effective masses. 

Equations (||)-(0) provide a connection between the total electron-hole density a-, and uncorrelated quasiparticle 
densities n^ and n^. The quasiparticle densities in turn govern the screeningl3 and therefore the strength of 
interaction between particles, which defines uniquely the set of binding energies and scattering phase shifts which 
enter Eqs. (|5|) and (|6|) for the two-body partition functions. These partition functions in turn define the ratio between 
n° and n^"'"'" via Eq. (Q). Thus, to find the degree of ionization of the electron- hole plasma one must solve the 
system of Eqs. @-(|3) self-consistently, together with a reasonable model of the screened interaction. 

In the next section we discuss the statically screened Coulomb potential which we use to model the interaction 
between-4)articles in an exciton/electron-hole plasma, and present results from the application of the variable phase 
methodE3 to scattering and bound states in this potential. In Section III wc present and discuss the results of 



calculations of partition functions and the degree of ionization of the electron-hole plasma. In Appendix H we 
derive the basic equations of the variable phase method, which is used for calculation of scattering phase shifts and 
binding energies. The 2D analogue of the Beth-Uhlenbeck formula and the modified law of mass action are derived 
in Appendix ^. 

II. STATICALLY SCREENED COULOMB POTENTIAL 

There is an extensive literature deahng with different aspects of the screened Coulomb interaction in 2D systems.i 
In this paper we model this interaction by the well-known Thomas- Fermi expression for a statically screened Coulomb 
potentialO 

= T^|^-|9.[Hofep)-A^ofep)]} , (7) 

where gs is the 2D screening wavenumber (which depends on temperature and carrier density), e is the static 
dielectric constant of the semiconductor, Jq(x)^ A^o(2;), and Ho(a;) are the Bessel, Neumann, and Struve functions 
respectively. The upper sign in Eq. (|^) is for electron-hole attraction, the lower sign is for electron-electron or hole-hole 
repulsion. nn 

Being the long-wavelength, static limit of the random phase approximation for a purely 2D case,c3E3 Eq. (|^) is 
the simplest model for the screened Coulomb potential in 2D. Nevertheless, this expression reflects the fact that 
the statically screened potential in 2D decreases at large distances ^loypf than in the 3D case (as a power law 
rather than exponentially). Despite numerous realistic corrections JH3c3'lZI Eq. (Hi jfemains the most widely used 
approximation for the 2D screening, especially for the screened exciton problem Su^lJ Optically active [ra = 0) 
bound states in the-attLactive, statically screened Coulomb potential [upper signjp Eq. (Q)] have been studied using a 
variational metliod,L3'Ej by a numerical pmcedure based on a shooting method]^ and more recently using the WKB 
approximatiorCj and perturbation theoryJlj As mentioned above, for the partition function calculation, all states are 
needed, bound and unbound, optically active and inactive. None of the above methods is suitable for analysis of 
shallow bound states and low-energy scattering states. 

We use fot| calculation of the scattering P|ba,se shifts and bound state energies entering Eqs. (0) and (pi) the 2D 
formulationuj of the variable phase method. E2l In this method the scattering phase shift and the function defining 
bound-state energies can be obtained as the large distance limit of the phase function, which satisfies the first-order, 
nonlinear Riccati equation originating from the radial Schrodinger equation (see Appendix H). The variable phase 
method is especially effective for calculation of the shallow-state binding energies and low-energy scattering phase 
shifts. 

In Fig. |l we show the /c-dependence of the scattering phase shifts for the attractive and repulsive Thomas-Fermi 
potentials [both signs of Eq. (1^)] with the screening wavenumber qs = 0.2/a*, where a* — eh? /{jie?) is the 3D 
exciton Bohr radius. The scattering phase shifts are negative for the repulsive potential and positive for the attractive 
potential. For the repulsive potential all zero-energy phase shifts vanish, 5m{k = Q) = for all angular momenta 
771. For the attractive potential 

hm 5„i{k) ^ v^TT , (8) 



where Vm is the number of bound states.Ej Equation (ph is the 2D analogue of Levinson's theoremL^Efl which 
connects the zero-energy scattering phase shift with the number of bound states for non-relativistic particles in 3D. 



This theorem has been kpojSHij^r almost five decades, however its applicabihty to the 2D scattering problem has been 
considered only recGBtly.E3cjE3 

We have recentlyt3 used Levinson's theorem in the form of Eq. (|^) for bound-state counting in the attractive 
Thomas- Fermi potential, Eq. ([7|) , and found a remarkably simple relation between the number of bound states and 
the screening wavenumber q^. Witk-decreasing screening, new bound states appear at critical values of the screening 
length given by the simple formulae 



1 



(2|m| + iy- 1)(2|to| + v) 



i^=l,2, 



(9) 



Equation (0) can be easily inverted and the number of bound states for given m and qs can be expressed as 

Vm = max{0 , vq - 2|m|} , (10) 

where 



vq = 



v/8/fea*) + 1 + 1 



(11) 



= 0. Here, and in Eqs. (^^ and (|14[), the bold square brackets designaie 
qs, Eq. (pH) gives a number 2.5 times smaller than the WKB ^atimate*^^' 



is the number of bound states with m = 

the integer part of a number. For small 

for the maximum number of bound s-states. The Bargmann bound conditiorO (re-stated for the 2D caseEj) for the 

attractive potential (H) is z/„j < 1 / {mqsa*) . This was also found to give a gross over-estimate of the number of bound 

states. 

The total number of bound states, Nb, for a given g^a* 



|m|„ 



Nt = vo + 2 Y^ v„ 



(12) 



can also be found explicitly as follows. From Eq. (H) the maximum possible value of \m\ for the state which remains 
bound is 



771 h 



v/8/(g,a*) + 1 -1 



l/Q -1 



(13) 



Then the sum in Eq. (n2|) can be easily evaluated using Eqs. (Kw and dl3): 

Nb = lyo + 2\m\„iax{l^Q - \m\max - 1) 



1^0+2 



Z/Q - 1 



f - 1 



vo-l 



(14) 



For small qsO,* (weakly screened potential) a simple approximate expression for the total number of bound states 
follows from substitution of Eq. (pTJ) into Eq. (Ill): 



Nb ^ -^ 



1 



Qsa* 



(15) 



Thus, for the weakly screened Thomas-Fermi potential the bound-state sum in the partition function, Eq. (g|), has 
a finite number of terms which is approximately equal to tlae. screening radius 1/qs measured in units of the Bohr 
radius. The WKB estimate of the number of bound states^ gives a different (square root) dependence of Nb on 
l/{qsa*) for small q^a*. The reason for this difference is that in Ref. p^ only m ~ states are considered, whereas 
all values of m are needed to obtain the result of Eq. ( p^ . 

As the screening is reduced Nb given by Eq. (O) exhibits steps of ever increasing height. In order for the limit 
of Eq. ( p^ ) to be meaningful the step height should be smaller than Nb itself, i.e., the normalized number of 
bound states, Nb/N^"^ ~ {qsa*)Nb, should converge to unity as g^a* — » 0. As can be seen in Fig. ^ this number 
oscillates around unity with the amplitude of oscillations decreasing with increasing 1/ [qsa*). It can be shown that 
for qsO* — > the amplitude of these oscillations is proportional to (gga*)^" and their period is proportional to 



In order to calculate the partition function, the bound state energies are required. T hese can also be obtained using 
the variable phase method, and the necessary equations are presented in Appendix [A 4 Numerical results for the 
attractive screened Coulomb potential [upper sign in Eq. (Q)] are presented in Fig. ||. In this figure the energies Em.u 
of the several lowest bound states of the screened exciton are shown as a function of the screening wavenumber gsO*. 
Here the energies are measured in effective exciton Rydberg [Ry* = Ti /(2/ie/ia* )] and we use the same classification 
of energy levels as in Ref. |3^, i.e., each energy level is characterized by the angular momentum quantum number ra 
and another number v which numerates different bound states for a given 771, with [v — 1) being the number of 
non-zero nodes of the radial wave function. For m = (s-states) the calculated energies are consistent with those 
obtained by J. Lee et oO using a variational method. 

III. PARTITION FUNCTIONS AND IONIZATION DEGREE 

Before we present the results of calculations of the partition functions and the ionization degree, we would like to 
discuss an important consequence of Levinson's theorem for the statistical mechanics of the 2D gas with an attractive 
interaction between its particles. The bound-state sum Zbound = X^m u 6xp(— /3i?m,i/) entering the two-body partition 
function in Eq. (0) exhibits jumps whenever bound states disappear with increasing screening. We will now show that 
these jumps do not give rise to unphysical discontinuities in the partition function if the scattering states are properly 
taken into account. 

Integrating by parts the scattering term and using Levinson's theorem in the form of Eq. (18) we can rewrite Eq. (1q) 
as 

Zeh =- ^ {exp(-/3S,„,,) - 1} 

+ -^ / I >^ <5™(fc)) exp(-fcV'7T)fc*, (16) 



/ ^ <5™(fc) exp(-fcV'7T) fc* , 

•^0 Vm=-oo / 



where q^ = 2^ehkBT /h . The modified bound-state sum [the first term in Eq. (ttq)] does not exhibit jumps whenever 
bound states disappear with increasing screening. For non-zero temperature the scattering integral [the second term 
in Eq. jig)] is also a smooth function of the interaction strength, which can be understood from Fig. ^. In this figure 
the scattering phase shift 82 is plotted as a function of k for several values of \/{qsa*) close to the critical value 
l/{qsa*) — 10 when the first bound state with m = 2 appears. One can see that although (5m (0) has a jump when 
qs passes a critical value this jump does not influence the value of the scattering integral if the thermal wavenumber 
qT is larger than the interval of k in which 5m{k) changes rapidly. As shown in Fig. H, when the bound state 
disappears the phase shift is affected only in an infinitesimally thin region around fc = 0. For any non-zero temperature 
this transition region makes no contribution to the phase-shift integral. Thus the electron-hole interaction part of the 
partition function given by Eg. (Hq) is a smooth function of the interaction strength as expected from the general 
thermoidjpamic argument E3o Similar cancellation of the bound state sum discontinuities for a 3D plasma is well 
known .cic3 

The results r^f— ithe calculation of the two-body interaction part of the partition function for the model 
semiconductor ,tj'Ell for which the assumption Me = Mh = l^^h is made, are presented in Fig. |^. Calculations 
are performed for two values of the ratio of fc^T to the bulk excitonic Rydberg, fc^T/Ry* = 1 (three upper curves) 
and fc^T/Ry* = 5 (three lower curves), which roughly correspond to ZnSe (or GaN) and GaAs at room tempera- 
ture. Solid lines show the bound-state sum, Zhound — J2m v ^^P{~E!m,v/kBT), which exhibits jumps whenever bound 
states disappear with increasing screening. The electron-hole part of the partition function, Z^h, which is shown by 
dashed lines, is a smooth function of the screening parameter, the bound state sum discontinuities are compensated 
by the scattering state contributions. Dot-dashed lines show the sum Zeh + Zee, which enters the modified law of 
mass action, Eq. (Q) (when simplified for the model semiconductor). Note that the cancellation of the Zeh term by 
the Zee term for fc^T/Ry* = 5 is stronger than for fc^T/Ry* — 1. This can be explained by the enhanced role of 
scattering states for the higher ratio of fcsT to the excitonic Rydberg. The lower absolute value of Zeh + Zee ensures 
that thermodynamic properties of the 2D electron-hole plasma in GaAs are much closer to the ideal gas behavior 
than those in the case of the wide-gap semiconductor BJ 

For the two-component electron-hole plasma the Thomas-Fermi 2D screening wavenumber entering Eq. (0) is given 
in the Boltzmann limit hyEB 



qsd = 



'^^^^ / , o^ . R-y 



fiekksT 



{nl + nl) = 4.^(A*Vn°a*^). (17) 



iNote tnat we use uncorrelatea quasiparticle densities n'^ anaj-a^ tor tne calculation ot tne screening wavenumoer, 
since Eq. ( |17| ) is derived for the non-interacting 2D plasmaBJEj We assume that the screening by excitons is much 
smaher then the free carrier screening when exciton and free carrier densities are of the same order. However, if one 
calculates the screening wavenumber using the difference between the total number of carriers and the number of 
bound carriers, unphysical jumps appear in the dependence of the screening wavenumber on total density as shallow 
bound states disappear with increasing density. Thus, it is natural to calculate g^ on the basis of uncorrelated 
density (n^ + n^), which is a part of the total density behaving as an ideal gas (see Appendix ^ and which is a 
smooth function of the total density. 

For the model semiconductor, rS ^ n^ ^ an, and Eq. dl7^) can be further simplified to 



Ry* ^2 
QsU* = St: a- — — na* . (18) 

kBi 

Equation ( |l8| ) shows clearly the connection between the dimensionless screening parameter g^a* and the two main 
dimensionless parameters characterizing the 2D electron-hole plasma, namely the dimensionless density na* and 
temperature /csT/Ry*. In addition the role of the degree of ionization, a , introduced by Eq. (||) becomes more 
transparent. The parameter a enters Eq. (nSf) explicitly, governing the screening wavenumber which determines the 
strength of the interaction between charged particles in the plasma. In turn the degree of ionization itself depends on 
qsa* through the partition functions Z^h and Z^e- For the model semiconductor the modified law of mass action, 
Eq. (Q), can be rewritten as 

^co.r^*2 ^ 4^ {nla*^) ' ^ {Z,n + ^ee) , (19) 

and using Eqs. (|^) and (|l^ we get the following expression for the degree of ionization a: 

a = |l + Ml(z,, + Z,,)| . (20) 

Equation (^0|) together with Eqs. (||) and (H) for the partition functions allows us to calculate the degree of ionization 
a of the dilute (nondegenerate) 2D electron- hole plasma as a function of the screening parameter g^a*. The connection 
between qgU* and the total electron (hole) density n [Eq. (p^)] can be used for self-consistent calculations of a 
as functions of n for different temperatures or material parameters. The result of these calculations are shown 
in Fig. 0. Calculations are performed for the model semiconductor with the exciton Bohr radius a* and effective 
Rydberg Ry* corresponding to ZnSe (Fig. ga) and GaAs (Fig. ||b) and for room temperature {ksT = 300 K). The 
arrows indicate the points of crossover from Boltzmann to Fermi statistics, n — 2/A^^ . As mentioned in Section |, 
the nondegenerate treatment is more adequate for the wide-gap material. 

On the same plot we show by the dashed lines the degree of ionization calculated using a simple law of mass 
action with a single bound state (the ground state of the screened exciton) . It can be seen from the figure that the 
degree of ionization is well described by the single-bound-state mass action law only for low carrier densitiesO For 
high densities (but remaining in the nondegenerate regime) the role of scattering states becomes essential. Instead 
of the unphysical behavior predicted by the simple mass action law, in which the degree of ionization decreases with 
increasing density, we find that the degree of ionization increases at higher densities. 

A minimum on the curve showing the density dependence of the degree of ionization has the following explanation. 
At low densities the main contribution to the correlated density comes from the ground exciton state, which is almost 
unscreened. This state in 2D is at least nine times deeper than the first excited state. Therefore, the simple single- 
bound-state law of mass action is a good approximation at low densities, but not as n ^ - when the number of 
bound states becomes larger than the ground state contribution to the partition function, see Ref. ^. The standard 
law of mass action states that the density of bound states is proportional to uncorrelated density squared, which 
reflects the fact that at fixed temperature (room temperature in our case) and at low density most of carriers occupy 
the high-energy ionized states in the continuum rather than the bound states. The low-density high-temperature 
electron-hole plasma behaves as an ideal gas with the degree of ionization close to unity. Thus, at low density 
the correlated density is proportional to the square of the total density and the degree of ionization decreases with 
increasing density. However, with the further increase in the total density screening becomes important and the inter- 
particle correlation caused by the Coulomb interaction starts to decrease. Correspondingly, the degree of ionization 
changes the character of its density dependence. There is a certain value of density, which corresponds to the minimal 
value of the degree of ionization. 

As expected, in wide-gap semiconductors the calculated degree of ionization is much lower than in GaAs for the 
same temperature and carrier density. For both materials the calculated degree of ionization of the room-temperature 



2D electron-hole plasma reaches its minimum at a certain density. The same happens for a 3D plasmajl3e3 however 
the minimal value of the degree of ionization for the 3D plasma is much higher than in the 2D case (compare Fig. g 
with Fig. 1 in Ref. Q). This is due to the much enhanced binding energy in 2D. 

The inclusion of Fermi statistics and phase-space filling, which is beyond the scope of the present paper, would 
provide a sharper rise of a at high carrier densities as the phase space available for the construction of exciton 
states is restricted. This will apply to both wide-gap and narrow-gap semiconductors. In the foregoing discussion 
we have assumed a purely 2D plasma. This assumption gives an overall overestimate of exciton binding energies, 
compared to a real, finite-width quantum well, forj-which unscreened exciton binding energies are lower and the finite 
thickness correction enhances the screening effectEZl Thus, the results shown in Fig. O should be considered as lower 
bound estimates of the degree of ionization of the electron-hole plasmas in ZnSe and GaAs quantum wells at room 
temperature. Note, that even this lower estimate does not give a value of the degree of ionization of a plasma in 
a ZnSe quantum well below 0.33. It means that at room temperature at least one third of the carriers are always 
unbound, which has to be taken into account in gain calculations. 

IV. CONCLUSION 

We have calculated the degree of ionization of the 2D electron-hole plasma taking into account all screened exciton 
bound states as well as scattering states. It has been shown that the scattering state contribution changes the 
character of the density dependence of the degree of ionization. We have found that the degree of ionization of the 
2D plasma reaches its minimal value at intermediate densities and approaches unity at high densities, which differs 
from the result based on the simple law of mass action. 

The calculated degree of ionization of the electron-hole plasma in a ZnSe quantum well is significantly lower than 
in a GaAs quantum well with the same carrier density and temperature. Therefore, excitonic processes should be 
considered for gain calculations in quantum wells based on wide-gap semiconductors. However, at room temperature 
at least one third of the carriers in ZnSe wells is shown to be unbound, which allows us to speculate that the most 
likely lasing mechanism at moderate density is exciton/free-carrier scattering. 

Most of the results presented here are obtained for the model system with equal electron and hole effective masses. 
For wide-gap semiconductors at room temperature {ksT ^ 1 Ry*) this approximation is good, since Z^e is much 
smaller than Z^i^ and the influence of the electron-electron part of the partition function on the degree of ionization 
is not significant. In the case of an extreme difference between electron and hole masses the model fails, e.g., lighter 
quasiparticles can be degenerate, when heavy quasiparticles are nondegenerate. 

The variable phase method is a powerful tool for studying scattering and bound states in any short-range potential. 
This method enabled us to find hitherto undiscovered properties of a Coulomb potential statically screened by a 
2D electron gas. The same approach can be applied to a more refined potential, which takes into account Friedel 
oscillations and the finite thickness of the 2D layer. 
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APPENDIX A: VARIABLE PHASE METHOD IN TWO DIMENSIONS 

In this Appendix we derive the basic equations of the variaMe phase approach in two dimensions from the radial 
Schrodinger equation. This derivation is similar to that in 3D.l3 

1. Scattering phase shifts 

The relative in-plane motion of two interacting particles with masses Ma and Mf, and the energy of relative 
motion E can be considered as a motion of a particle with the mass fXat — MaMh/{Ma + Mb) and energy E, 
moving in an external central potential V{p). This motion is described by the wave function satisfying the stationary 
Schrodinger equation 



Hreli^ 
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i}j + v{p)i: = ^v 



(Al) 



2^ab \pdp dp P^ dip^ ^ 

Owing to the axial symmetry of the potential V{p), we can separate variables in the expression for the wave function 

V'™(p,</?) = i?™(p)e™'^, m = 0,±l,±2, ... . (A2) 

The equation for the radial function Rm{p) reads 



R„ 



-R'^ +{e- U{p) -^]Rm = 



(A3) 



where k^ = 2pabE/h^ and U{p) — 2p,abV{p)/h^ . In what follows we consider ?ti > only, as R^rn{p) — Rm{p)- 

We assume that the interaction potential vanishes at infinity (the precise decay rate will be discussed later). Then 
at large distances the radial function satisfies the free Bessel equation, whose general solution is 



Rm (p) = Am [J„i (kp) cos S,n ~ N,n (kp) siu S„ 

1/2 



p— *oo 



A, 



2 
nkp 



cos{kp - (2m + l)7r/4 + S^) 



(A4) 



where Sm is the scattering phase shift,EZrE3 Jm{kp) and Nm{kp) are the Bessel and the Neumann functions, 
respectively. 

In the variable phase approach. Am and Sm are considered not as constants but as functions of the distance p. 
The amplitude function Am{p) and the phase function Sm{p) are introduced by the equation 



Rm{p) = Am{p)[Jm{kp) COS Smip) - Nm{kp)sin6mip)] 

with the additional condition, which we are free to choose as 

R'm (P) = Am {p) [J'm i^p) COS 5m {p) - N^ (kp) siu 5m (p)] 



(A5) 



(A6) 



where the prime indicates differentiation with respect to p. The phase function 5m{p) has a natural physical 

interpretation as being the phase shift produced by a potential cut-off at a distance p. 

Differentiating Eq. ( |A6| ) and substituting the resulting expression, together with Eqs. ( [A^ ) and (A6), into Eq. ( [A^ ) 
we get 



A'm{p)[JLikp)cos5mip) - Nl^{kp)smSm{p)] 
-5',Jp)Amip)[Jmikp) sin5m{p) + N'^{kp) cos5m{p)] 

= U{p)Am{p)[Jm{kp) COS 5m{p) " N m{k p) Sin 5 m{p)] 



(A7) 



To obtain Eq. (A7) we used that the functions Jm{kp) and Nm{kp) satisfy the free Bessel equation: 



F" +-F' 

p 



k-—\ Fm = 



Equating the derivative of Eq. ( [A5| ) to Eq. ( A6 ) implies the following condition on the derivatives of the amplitude 
and the phase functions: 



A'm (P) [Jm (kp) cos 5m (p) - Nm {kp) siu 5m (p)] 

= 5'^{p)Am{p)[J7n{kp) sinSmip) + Nm{kp) COS 5m{p)] 



Substituting A'^{p), obtained from Eq. ([AS|), in Eq. ( [AT] ) yields 

- 5'm{p)[Jm{kp)NUkp) - Nm{kp)j'^{kp)] 

= U{p)[Jm{kp) COS 5m{p) - Nm{kp) sin 5m{p)f ■ 
Equation (Ra) can be simplified further, using the Wronskian of the Bessel functions 

W{Jm{x),Nm{x)} = Jm{x)-r-Nm{x) ~ Nm{x)-—Jm{x) = 

ax ax TTx 



(A8) 



(A9) 



and thus becomes 



dp 



X [Jm{kp) COS Sjn{p) ~ ^m(fcp) Slli Smip)Y 



(AlO) 



This phase equation, Eq. (AlC), is a first-order, non- linear differential equation of the Ricatti type, which must be 
solved with the initial condition 



<5m(0) 







(All) 



thus ensuring that the radial function does not diverge at p = 0. The total scattering phase shift Sm can be obtained 
as a large distance limit of the phase function Sm (p) : 



lim Smip) 



p — >oo 



For numerical convenience, instead of the initial condition Eq. ( All ), the small-p expansion is used 



6mip) 



22'"+i(m!)2 



TT^ /\ /2m+l , I 

U{p)p dp , p 







(A12) 



(A13) 



From Eq. ( |A10| ) and the asymptotic expansions of the Bessel functions one can see that the variable phase method 
is applicable only if the scattering potential U{p) satisfies the necessary conditions: 



U{p') dp' -^ 



oo , 



(A14) 



and 



p'C/(p) 



(A15) 



The statically screened Coulomb potential Vs (p) , defined by Eq. (^ , behaves like p ^ at small distances and like p '^ 
at large distances. Such behavior allows the application of the variable phase method to this potential. 

2. Bound state energies 

For the states with negative energy of the relative motion (bound states), the wavenumber k is imaginary, k = in, 
and we introduce the function rj^ip, k) vanishing in the origin and satisfying a non-linear equation 



^pVn.{p,n) 



^pu{p) 



ImiKp) cos77„(p, k) + -Krn{np) sin?7,„(p, k) 



(A16) 



where Imi^p) and Km{np) arc the modifie d Bes sel functions of the first and second kind, respectively. Equation 
( [Al(j ) is derived in the same fashion as Eq. ( AlO ). The fmictions I„i{np) and Kminp) represent two linearly 
independent solutions of the free radial- wave Schrodingcr equation for the negative value of energy, E — —fi k^ I2p,ab, 
and cot 77m characterizes the weights of the diverging [/^(Kp)] and converging [Km{np)\ solutions as p — > cx) . For 
the bound state, the diverging solution vanishes, implying the asymptotic condition 



r]rn{p ^ oo, K^) = {v-1/2)tt, 1/^1,2, 



(A17) 



Here v numerates the bound states for a giv en m and (z^ — 1) is the number of non-zero nodes of the radial wave 
function. For numerical solution of Eq. (AlC), instead of the boundary condition 77™ (0, k) — 0, an asymptotic initial 
condition [analogous to the condition Eq. (A13) for the phase function Sm{p)] is used. 



APPENDIX B: BETH-UHLENBECK FORMULA IN TWO DIMENSIONS 



In this Appendix we derive Eq. (|h, which is the 2D analogue of the Beth-Uhlenbeck formula,li3 and the modified 
law of the mass action, Eq. (Q). This derivation is similar to the|-analysis used for the calculation of the second virial 
coefficient of low-density '^He and *He monolayers on graphite.L3 

Let us consider a binary mixture of components a and h in two dimensions. The grand partition function of the 
mixture is given by 



n{za, zb, A,T) ^ Y. Qn^,n, z^- zf " , 



(Bl) 



Na.,Ni 



where Za and Zb are the fugacities {za = e'^^" , with /ia being the chemical potential of the component a), A is 
the area of the 2D system, and QjVa.JVi, is the partition function defined as 



iN^,Ni 



{A,T) = Tre" 



-0H{N^.Nt) 



(B2) 



where the trace is to be taken over all states of the system that has Na particles of the type a and Ni, particles 
of the type b in the area A. 

We now expand the quantity InQ as a power series in Zq, Zb 



la,h 



\nn{za,Zb,A,T) ^ Y, ACi^,i, z[- zl' , 
The density of the component a is given by 



1 '^ 1 o 

— Za 3— lnS2 

A dZa 



la lb 



(B3) 



(B4) 



la,lb 



From this point we consider the low-density limit, Za, Zb <SC 1, and neglect all the terms higher than z^ in Eq. 
Then 

ria ~ Ci.QZa + Ci^iZaZb + 2C2fizl . 



From comparing corresponding powers in Eqs. (B3) and (Bl) we get 

Ci,o - Qi.oM , 



(B5) 
(B6) 



Cia = (Qia-Qi,oQo,i)M: 



and 



^2,0 



Q2,o-2 2lo)M 



(B7) 



(B8) 



The next step is to calculate the partition functions entering Eqs. (|B6|)-(P38|). First of all the one-particle partition 
function, Qi,0) is given by 



Q.o(A,T).,.^/d^kexp(-^^) 



(2^)2 J - V '"2M„, 

where Qa is a quantum state degeneracy and Xm^ is a thermal wavelength, 

2-K(3h^ 



A 



\^ — 



M„. 



This yields 



^hO - T2- 



(B9) 



(BIO) 



(Bll) 



'Ma 



10 



In order to find the two-particle partition function, Qi,i, it is useful to separate the center of mass motion and the 
relative motion of the two particles: 



■^Ma+Mb 



A 



gagt 



X2 

'^Ma + Mb 



Jd'pJ2\Mp)\^e- 



■I3E„ 



(B12) 



where the factor A/\\j , ^^ appears from performing the summation over all center-of-mass momenta, the Hamil- 

tonian Hrei of the relative motion is given in Appendix ^ and the sum in Eq. (B12) is taken over all different 
solutions of Eq. ( Al ) . 



For the corresponding two-body system of noninteracting distinguishable particles, one would have 



(B13) 



where the superscript ^^' refers to quantities of the noninteracting system. The two-body interaction part of the 
partition function is then defined by 



-' n 



(B14) 



Thus 



Qi.i(AT) = Qf}{A,T) + gagt-^^^ Zab 



A 



(B15) 



Ma+Mb 



r(0) 



To analyze Eq. ( B14 ) further we must study the energy spectra En and En- For the noninteracting system. 
En forms a continuum. We write 



Ei'^ = 



2 1,2 






(B16) 



which defines the relative wave number k. Then for the system of two noninteracting distinguishable particles the 
function Q[ { given by Eq. ( B13 ) can be easily evaluated as 



Qf}{A,T) ^ gagb^ 



A 



A A 



d k exp -/3 



2/-iafc 



gagb 



X2 \2 

^M^ + Mb ^tJ.ab 



ga 



A 

12" 



gb 



A 

^Mb 



Qi,oQo,- 



(B17) 



For the interacting system, the spectrum of En in general contains a discrete set of values Eb, corresponding to 
two-body bound states, and a continuum. In the continuum, we define the wave number k for the interacting system 
by putting 



En - 



2Hab 



(B18) 



Let g(k) dk be the number of states with wave number lying betwe en k and k + dk , and let g^'^'{k) dk denote 
the corresponding quantity for the noninteracting system. Then Eq. (B14) can be written in the form 



Zab = Y.'rf'^^+ dk{g(k)-g^°\k)]e^^(-(3^ 

B Jo \ ^Ma 



2 (.2 
b 



(B19) 
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The difference in density of states is related to tfie scattering pfiase shifts by the following argument. cS The relative 
wave function can be factorized [see Eg. (|A2| )] into a product of a trivial azimuthal part and non-trivial radial wave 
function Rm{p), which satisfies Eq. (E^). For large value of p where the potential is assumed negligible, 



Rm{p — > oo) oc cos{kp — (2m + l)7r/4 + Sm{k)} , 



(B20) 



which defines the phase shift Sm{k) of the m-th partial wave. For the noninteracting system all the phase shifts 
dm{k) = 0. If the system is placed within a circle of radius R, the vanishing of the wave function at the boundary 
requires that the allowed values of k are given by 



kR - (2m+l)7r/4 + Smik) 



for the interacting system, and 



kR - (2TO+l)7r/4 



n+ - \ TT 



71" , 



(B21) 



(B22) 



for the noninteracting system, where n = 0, 1, 2, , 
by the respective amounts Afc and AA;'*'^: 



A/c 



For a given m, changing n by one unit causes k to change 



TT 



R + {d5m{k)/dk) ' 



A/fc(o) = ^ 
R 



(B23) 
(B24) 



These are the spacings of eigenvalues for a given m. Let the number of states of a given m with wave number lying 
between k and k + dk be denoted by g,n{k)dk and gfn {k) dk for the two cases. We must have 

(B25) 
(B26) 



(B27) 
(B28) 







9m{k) 
9t\k) 


A/c = 
A/c = 


1, 
1 


9m{k) 


1 

TT 


R + 
V 


dSm {k) 
dk 



9l°Hk) 



R 



Therefore 



9m{k) - 9lil\k) = - 



1 ddm{k) 



TT dk 



Summing Eq. ( B29| ) over all allowed m we obtain 



Substituting Eq. ( ^ ) into Eq. ( |B19| ) yields 



Zeh = y^ ' 



-f}EB 



B 
1 

71" ./O 



f; *#))e.pr-„?!^u* 



V,?n— — oo 



dk 



"^Pab 



(B29) 



(B30) 



(B31) 



which coincides with Eq. (1q) if we change the notation for the bound-state energy from Eb to Em,i^, where subscript 
h' enumerates bound states with a given m. 

Now, having evaluated the interaction pa rt of the partition function, we can obtain the co efficie nt Ci,i needed in 
the density expansion, by substituting Eq. (B15) into Eq. (B7) and taking into account Eq. (B17), 
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Cl,l — gagb-r2 Zab — -^p^t, T2 — T2 — ^ab- (B32) 

Up to this point we have considered a system of two distinguishable particles (e.g. an electron and a hole). To 



calculate the remaining coefficient, C2,o, in Eq. (B5) we must now consider the interaction of indistinguishable 
particles. For such particles the sum over different m is modified, e.g., for two spinless bosons the relative wave 
function must be symmetric, this means that only even values of m are possible. For the case of two fermions with 
the same spin the coordinate wave function must by antisymmetric. Let us consider a case of the two fermions or 
bosons of the same type a , with the quantum state degeneracy ga- This degeneracy is usually associated with the 
particle spin s, so that ga = 2s + 1. For a two-particle system there are g^ = (2s + 1)^ spin states of which a 
fraction, (s + l)/(2s + 1) = {ga + l)/(25a)j are symmetric and s/(2s + 1) = {ga — l)/(2ga) are antisymmetric. For 
a fermion system, the symmetric spin states must be multiplied by antisymmetric spatial states (odd m) and the 
antisymmetric spin states multiply symmetric spatial states (even m). A similar argument is applicable for bosons. 
Thus, for the 2D system of interacting particles of the type a, the two-body partition function is 



A 
X 



Q2,o{AT) = Q^^l{A,T) + gl^Zaa , (B33) 



A /a 



where Zaa is defined by 



Zaa - 2^ 



ga ± {-ly 
ga 



g-/3£„ 



exp -13-—- dk , (B34) 



- E 



m— — oo 



ga Jo dk \ M, 



Here and in what follows the upper sign is for bosons and the lo wer sign is for fermions. For the case of repulsive 
fermions (which do not form bound states) and for g^ = 2 , Eq. ( B34 ) reduces to Eq. (||). 



Let us now calculate the two-body partition function of nonintcracting bosons or fermions, Qg o^ which is given by 



S^°^(AT) = T^f fd'pY^'{{ga±l)\4"KKp)\ 



^2Ma 

+ {ga T 1) l^r(k,p)p} exp (-/5^) dk , (B35) 

where ips (k, p) and ipa (k, p) are symmetric and antisymmetric eigenfunctions of the nonintcracting Hamiltonian 
of the relative motion, i.e., 

V'f(k,p) - yicos(k.p) (B36) 

and 

V'i°)(k,p) - y|sin(k-p). (B37) 



The summation in Eq. (B35) must be performed over all different two-particle states. Wave vectors k and — k 
correspond to the same state (k and — k transfer into each other upon exchange of the two indistinguishable 
particles), this state should not be counted twice and 



Y.' 



1 A 



2 (27r) 



d^k . (B38) 



Then Qj o <^an be written as 



Q^°l{A, T) - ga ^ [{ga ± 1)X, + (g, T l)^a} , (B39) 
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where 



1 /■ f fi k'' 

Is = j:^ / rf2prf^kcos2(k-p)expf -/3-^ 

1 / A 



2 V2A 



(B40) 



and 



la = 



i^ I dVrf^k sm2(k . p) cxp (^-/3^ 



(2 

1 / A 



2 X'^^M^ 



(B41) 



Finally, 



i(o) 



Q^:^(^,T) 



2 r Ai, 



1 A 

4 ^M„ 



r,2l,0 ± A 9a ^2 



1 

4 5" A 



A 



Afa 



(B42) 



Substituting Eqs. ( p33| ) and ( |B42D into Eq. @ yields 



^2,0 









4Ak 



(B43) 



Substituting Eqs. (Bll), (B32), and (B43) into Eq. (B5) we find the following expression for the total density of 
the component a: 



9a 



± I 3a 



2 ^Ma 



+ A 



9a 



.zl + A^ 



9a 9b 



tiab \2 



X2 

-^A/t, 



Zab Za Zb 



(B44) 



The first two terms in Eq. (B44) do not depend on interparticle interaction (however the second term depends via 
its sign on the statistics of the particles). It is natural to call the sum of two first terms in Eq. (B44) as uncorrelated 
density, we denote it raa n^. Note that these terms are the first two terms in the low-density expansion of the 
well-known expressiont^l for the density of the noninteracting 2D Bose or Fermi gases: 

TTl^ln(lT^a). 



A 



Afa 



The sum of the last two terms in Eq. (B44) is the interaction-dependent correlated density n™''''. 

In the low-de nsity limit the second term in Eq. (B44) is much smaller than the first one, and within the same 
accuracy as Eq. ( B44 ) we can write 



^nKA^.^ab 



(B45) 



Equation (B45) coincides with Eq. M), and it constitutes the modified law of the mass action in two dimensions. 
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FIG. 1. Scattering phase shifts versus the in-plane wave vector k (in units of inverse Bohr radius 1/a*) for a 2D particle 
in a screened Coulomb potential, Eq.(M). Screening wavenumber qs — Q.2/a* . For the attractive potential all phase shifts are 
positive and for the repulsive potential they are negative. Numbers show m values. 
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FIG. 2. The normalized number of bound states {qsa*)Nb as a function of the inverse screening parameter, l/{qsa*), for 
qsO.* < 1- 
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FIG. 3. The bound-state energies Em,i^ of the 2D exciton in exciton Rydberg units are shown as a function of the screening 
parameter qsU* for different m values. Solid lines show m — states (-Eo.i, £^0,2, and £0,3); dashed lines show m — 1 
states (i5i,i and Si, 2); the dot-dashed line shows the lowest state with m = 2 (£2,1). Inset is a blowup near g^a* — 0. 
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FIG. 4. The scattering phase shift & is shown as a function of the in-plane wave vector k (measured in inverse exciton Bohr 
radii) for several values of the inverse screening parameter close to the critical value, l/{qsa*) — 10. Solid line: l/{qsa*) = 9.9; 
dashed line: l/{qsa*) — 9.95; dot-dashed line l/{qsa) = 9.98. Dots show 52{k) for l/(gsa*) — 10.1 (a shallow bound state 
with m = 2 has just appeared). 
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FIG. 5. The two-body interaction part of the partition function versus the inverse screening parameter l/{qsa*) for two 
values of fcsT/Ry*. Three upper curves: ksT = 1 Ry*; three lower curves: fesT = 5 Ry*. Solid hues show the bound state 
contributions Zbound only; dashed lines show Z^h; dot-dashed lines show Zeh + Zee. 
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FIG. 6. The degree of ionization (solid lines) of the nondegenerate 2D electron-hole plasma as a function of the total electron 
density at room temperature, calculated for the model semiconductor with the effective Bohr radius and excitonic Rydberg of 
(a) ZnSe, and (b) GaAs. The arrows indicate n — 2/\\i^ for ZnSe and GaAs at room temperature. Dashed lines show the 
degree of ionization calculated using a simple law of mass action with a single bound state. 
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